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ified the Transfer Integral method to adapt it to the calculation of both extensive 
(partition function, entropy, specific heat, etc) and non-extensive (order parame- 
ter and correlation length) thermodynamic quantities of finite sequences with open 
boundary conditions, and applied the modified procedure to two different dynamical 
models. We showed that rounding of the transition clearly takes place when the 
length of the sequence is decreased. We also performed a finite-size scaling analysis 
of the two models and showed that the singular part of the free energy can indeed 
be expressed in terms of an homogeneous function. However, both the correlation 
length £ and the average separation between paired bases (y) diverge at the melting 
transition, so that it is no longer clear to which of these two quantities the length 
L of the system should be compared. Moreover, Josephson's identity is satisfied for 
none of the investigated models, so that the derivation of the characteristic expo- 
nents which appear, for example, in the expression of the specific heat, requires some 
care. 
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I. INTRODUCTION 



Real systems manifesting critical behavior have necessarily finite volume. However it is 
well-known that the finiteness of the system size lets the critical singularities disappear and 
smears out the phase transition (see for example [1-5]). The obvious argument which may 
reconcile these two aspects is that, as the finite size of the system is increased and passes 
through a critical value which characterizes the border of the thermodynamical domain, 
the thermodynamical limit is approximately reached and consequently the critical singu- 
larities manifest themselves. Thus it is crucial to thoroughly understand the evolution of 
singularities with respect to the system size and estimate the critical size above which the 
thermodynamical limit is attained. 

An efficient tool to analyze the volume dependence of a critical phenomenon is the finite 
size scaling theory. Besides providing information on the rounding of critical singularities 
and the shift in the critical point, this theory is also an alternative way to determine critical 
exponents characterizing the phase transition. 

Finite size scaling theory was developed by Fisher and Barber in the early seventies 
During the last thirty years it has been applied to various systems exhibiting both first 
and second order phase transitions. Among hundreds of subjects we can mention the study 
of finite size effects at first order transitions by gaussian approximation [7-8], the Gibbs 
ensemble 9|, five dimensional Ising model [10-11], percolation models [12-13], stochastic 
sandpiles 14j, six-dimensional Ising system 15[], Baxter- Wu model 16], two dimensional 
anisotropic Heisenberg model 171] ... 

In nature, the majority of phase transitions are sharp and discontinuous first order tran- 
sitions. Experimental UV absorption spectra of diluted DNA solutions reveal that the DNA 
melting transition belongs to this class. A widely used dynamical non-linear DNA model 
was proposed twenty years ago by Peyrard and Bishop This model involving harmonic 
interaction terms between successive base pairs was later improved by the contribution of 
Dauxois [l9] and the new model (DPB model) yields a sharp transition. We recently pro- 
posed an alternative DNA model (JB model) which respects the finiteness of the stacking 
energy, and showed that this model also exhibits a sharp first order phase transition 20]. 



Then in 



2l| we showed that for both models the generalized homogeneity assumption is not 



respected so that Josephson's identity (also known as the hyperscaling relation) is not valid. 
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We tentatively explained this fact by the divergence of the order parameter at the critical 
point. 

The goal of this article is to investigate the sequence length dependence of the DNA melt- 
ing transition. This is an important point since experiments dealing with DNA molecules 
are carried out with various sequence lengths. To this end we employ a modified transfer 
integral method adapted to finite chains with open boundary conditions. The two hamilto- 
nian DNA models to be studied, i.e. the DPB and the JB models, are briefly described in 
Sec. II. Section III deals with the transfer matrix theory for finite linear chains and Sec. IV 
is devoted to the finite size scaling analysis which leads to a better understanding of finite 
size effects. 



II. NON-LINEAR HAMILTONIAN MODELS FOR DNA 

The Hamiltonians of the two DNA models whose critical behaviour is studied in this 
article are of the form 

N ( 2 } 

H = + V M(Vn) + W(y n , y n „ x ) (1) 

n=l ^ ' 

where y n is the transverse stretching of the hydrogen bond between the nth pair of bases, 
while the one-particle Morse potential term 



V M (y n ) = D (1 - e'^f (2) 

models the binding energy of the same hydrogen bond. The choice of the nearest-neighbor 
interaction potential W(y n ,y n -i) is crucial since the type of the transition, which is a collec- 
tive effect, depends primarily on its form. The DPB model 19) assumes that the stacking 
energy is of the form 

W(y n ,y n ^) = j(y n -y n -i) 2 [l + pe -^+^] . (3) 

This non-linear stacking interaction has the particularity of having a coupling constant which 
drops from K(l + p) to K as the critical point is approached. This decreases the rigidity of 
the DNA chain close to the dissociation and yields a sharp, first-order transition. 
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This interaction potential still has the inconvenience that the stacking energy diverges 
when two paired bases separate. Taking into account the finiteness of the interaction between 
adjacent bases, we proposed a potential of the form [20] ] 

A H r i 

W(y n , y n - x ) = — [l - e - fe ("-) 2 J + K b (y n - y n _ x f (4) 

which, contrary to the model ([3]), depends only on the distance between base pairs. The 
small harmonic term, whose constant K b is 2000 times smaller than the parameter K of the 
DPB model, was introduced in order to take into account the stiffness of the phosphate-sugar 
backbone. 



Numerical values of the parameters are those of Refs. [19,20], that is D = 0.03 eV, 
a = 4.5A" 1 , a = 0.35A" 1 , K = O.OGeV A" 2 , p = 1 for the DPB model, and D = OMeV, 
a = 4.45 A -1 , AH = O.UeV, K b = 10- 5 eVA" 2 and b = 0.10 A' 2 for the JB model. 



III. TRANSFER INTEGRAL METHOD FOR FINITE CHAINS WITH OPEN 

BOUNDARY CONDITIONS 

A. The partition function and extensive thermodynamic quantitites 

Let us define the Transfer Integral (TI) kernel according to 

K(y n , y n ^) = exp [-/3 {V M {y n )/2 + V M {y n -i)/2 + W{y n , y n ^)}} (5) 

where (3 = {ksT)" 1 is the inverse temperature. In the following analysis we deal only with 
sequences having open boundary conditions. Then the partition function of the system can 
be expressed as 

Z = J d yi dy 2 ..- dy N e-^^l 2 K (y 2 , y x ) K (y 3 , y 2 ) ■ ■ ■ K (y N , y N ^) e^^ 2 . (6) 
The TI method consists in expanding the kernel of Eq. ([5]) in an orthonormal basis 

K (y n , = KMyn)MVn-i) (7) 

i 

where the and {A,} are the eigenvalues and eigenvectors of the integral operator and 
satisfy 
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J dx K(x, y) $i(x) = \<bi{y) . (8) 

This integral equation was solved by diagonalizing the symmetric TI operator K (x, y) on a 
regularly spaced grid defined between y min = — 200/a and y ma x = 4000/a with 1/a intervals. 
Numerical integrations were performed on the same grid. 

In this study we extended the transfer matrix approach for open chains developed in 22 j 



to adapt it to the calculation of the order parameter (y) and the correlation length £. Let 
us first consider extensive thermodynamic quantities. By introducing 

dye-^^Qiiy) (9) 
and by substituting the kernel expansion of Eq. ([7]) into Eq. ([6]), we get 

z = J2^' 1 - (10) 

i 

Determination of the partition function then allows the computation of extensive quantities 
of the system such as the free energy, the entropy and the specific heat : 



F = -k B T1n(Z) 
dF 
df 



dF 

S = -7^ (") 



d 2 F 

°y - ~ T dfi ■ 

In the thermodynamical limit N — > oo the major contribution to the partition function 
arises from the largest eigenvalue Ai and in this limit it is reasonable to drop the eigenvalues 
with % > 2. Neverthless, we will consider large DNA molecules as well as small ones. 
Consequently as many eigenvalues as possible must be taken into account in numerical 
computations. From the practical point of view, it was found that considering the first 400 
eigenvalues is enough to insure numerical convergence of the results presented below. 

B. The order parameter and the correlation length 



The order parameter of DNA melting transition is the mean separation of the bases 
averaged over the sites of the sequence : 
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1 - 

(y) = ? EW- ^ 



N 

71=1 



In order to reduce (?/„) to a form depending only on the eigenvalues and eigenvectors of the 
TI operator, we first write it as 

(Vn) = ^Jd yi dy 2 --- dy N y n e-e y ^ 2 K (y 2 , Vl ) K (y 3 , Jfe) • • • K {y N) y N - X ) e^ v ^' 2 . 

(13) 

Substituting Eq. (J7J) into Eq. (TI3~1) and defining 



6, = y dye-^MMftQ^y 
Y^ = J dy^(y)y^(y), 



(14) 



we obtain 



Gft> = W^^a^Af- 1 (15) 



and 



(Vn) = ^J2a i a j Y^Xr 1 \t n ( 16 ) 

U 

for n ^ 1,N. By evaluating the geometric summation that appears in Eq. flT2l) when 
replacing the (y n ) by their expressions in Eqs. (|T5|) - (|T6|) . we finally get 

2 N 

with rjj = Xi/Xj. The computation of the correlation length £ proceeds along similar lines 
although it is more elaborate. A sketch of the derivation and the analytical result can be 
found in Appendix A. 
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IV. FINITE-SIZE EFFECTS NEAR THE CRITICAL POINT 
A. Rounding of the melting transition of DNA 

It is well known that a finite-size system does not exhibit any phase transition. At the 
critical point its free energy is analytic and consequently all thermodynamical quantities are 
regular. Let L be the size of a system having a critical behaviour in the thermodynamical 
limit L — > oo. For this system finite-size effects manifest themselves as e~ L ^, where £ is the 
correlation length, by rounding the critical point singularity. In other words they become 
important over a region for which £ ~ L. A simple example of this rounding phenomena 
concerning the Ising model can be found in For an infinite-size Ising system, as the 
magnetic field H varies, the order parameter jumps discontinuously from — M cr to +M cr at 
the critical point H = 0. On the other hand if the system's size is finite, this transition occurs 
on a finite region of order AH ~ k^T ' j '(M cr L d ) with a large but finite slope ~ M\L d /(fcgT) 
where Ml is the most probable value of the magnetization in the finite system. 

For the two DNA models sketched in Sect. 2, the size L of the system is equal to the 
number N of base pairs in the sequence times the distance between two successive base 
pairs. This latter quantity playing no role in the dynamics of the investigated models, we 
will henceforth use indistinctly iV or L to refer to the size of the sequence. 

Given a sequence of length N, the first task consists in determining its critical tempera- 
ture, which we denote by T C (N). Among the several methods listed for example in (25j, we 
found it rather simple and convenient to search for the maximum of the specific heat Cy, 
which is more pronounced than that of the correlation length £, thus allowing for a more 
accurate localization of the temperature. Two observations confirm a posteriori that the 
critical temperatures thus obtained are correct. First, the shift in critical temperature is 
found to vary as a power of N, as predicted by finite-size scaling theory (see below). The 
top plot of Fig. 1 shows for example that T c — T C (N), where T c stands for T c (oo), varies as 
iV~ 1,00 and A -1,05 for the DPB and JB models, respectively. Note that this scaling is also in 
excellent agreement with the semi-empirical formula used by experimentalists to calculate 
the melting temperature of finite sequences. Moreover, as will be seen later (Figs. 3 and 
4) the curves for the temperature evolution of Cy, £, (y), etc... for sequences with differ- 
ent lengths N all coincide sufficiently far from the critical temperature when plotted as a 
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function of the reduced temperature 

, Art T-T C (N) 

In order to illustrate finite-size effects acting on DNA melting, we first computed the entropy 
per base, s = S/N, for an infinite chain and a short DNA sequence for both the DPB and 
the JB models. Results are shown in Fig. 2. At the thermodynamic limit, the entropy 
s is clearly discontinuous at the critical temperature, as is expected for first order phase 
transitions. In contrast, smooth curves are observed over the whole temperature range for 
the sequence with N = 100. We next computed the specific heat per base, cy = Cy/N, 
for increasing sequence length and temperature. The top and bottom plots of Fig. 3 show 
the temperature evolution of cy for seven values of N ranging from 100 to infinity for the 
DPB and JB models, respectively. It is seen in this figure that rounding manifests itself 
through a decrease in the maximum of cy as N decreases, but also through the fact that 
the sharp rise of cy takes place further and further from the critical temperature, that is, 
at increasingly larger values of |£(iV)|. This is particularly clear for the DPB model, which 
at the thermodynamic limit undergoes a very sharp transition, i.e. a transition that is 



noticeable only at very small values of \t\ = \t(oo)\ 2lj. Quite interestingly, examination of 



Fig. 3 also indicates that the two models consequently give very comparable results up to 
iV ~ 1000, while the narrower nature of the phase transition for the DPB model becomes 
apparent for longer sequences. 

At this point, it should be emphasized that boundary effects may become important 
when the size of the system is small. In order to check whether such boundary effects play a 
role in the results presented above, we repeated these computations with periodic boundary 
conditions instead of open ones and found that this alters only very little the results for s 
and cy down to N = 100. Conclusion is therefore that boundary effects play only a marginal 
role down to this size. 

Finally, we computed, for the JB model, the temperature evolution of the correlation 
length £ (Eq. (1A2I) ) and the order parameter (y) (Eq. (flTl) ) for increasing values of N. 
As the critical temperature is approached, the correlation length £ of a finite-size system 
is expected to increase according to the power law £ oc t u (N) till it reaches the system's 
dimension L and freezes. This behaviour can be checked in the middle plot of Fig. 1, which 
shows the evolution of the maximum of the correlation length (in units of the separation 
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between successive base pairs) as a function of N : the maximum of £ is indeed of the same 
order of magnitude as N and the curve scales as iV a97 . An exception however occurs for the 
last three points with N > 3000 : we will come back later to this point. The bottom plot of 
Fig. 4 additionally shows the temperature evolution of £ for seven values of iV ranging from 
100 to infinity. One observes just the same rounding effects as for the specific heat in Fig. 
3. This is again the case for the temperature evolution of the order parameter (y), which 
is drawn in the top plot of Fig. 4. This latter plot however displays a remarkable feature, 
in the sense that all the curves converge to the same limit at T C (N). To understand why 
this is the case, it must be realized that (y), the average separation between paired bases, 
is the only quantity which diverges at the critical temperature whatever the size N of the 
sequence, while, for example, cy and £ diverge for infinitely long chains but remain finite 
for finite chains. The limit towards which all curves converge in the top plot of Fig. 4 is 
thus just the approximation of infinity imposed by the numerical procedure (size of the grid, 
etc.). 



B. Finite-size scaling analysis 

The basic idea of finite-size scaling is that the correlation length £ is the only length 
that matters close to the critical temperature and that one just needs to compare the linear 
dimension L of the system to £ : rounding and shifting indeed set in as soon as L/£ ~ 1 . 
Since, by definition of the critical exponent v, £ grows as t~ u , one has (L/£) oc (tL 1 ^ u ) v as 
L — > oo or, L being proportional to the number N of paired bases, (L/£) oc (tN l / v ) v . In the 
absence of external field, it is therefore natural to write the singular part of the free energy 
of the finite-size system in the form 

fsin g = N- d Y{tN l ' v ). (18) 

where Y is some homogeneous function. Differentiating Eq. (ji~8~l) twice with respect to t, 
one obtains that cy is equal to 

c v = N p G(tN a ) , (19) 

where 
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2 , 
p = a 



1 

a = — 
v 



(20) 



and G is an homogeneous function which is proportional to the second derivative of Y. By 
using Josephson's identity, 2 — a = ud, where a is the critical exponent for cy (cy oc t~ a ), 
coefficients p and a can be recast in the form 



a 

P 

(21) 



v 
1 



cr = — 

v 



Conversely, if there occur several lengths that diverge at the critical point, as is the case for 
DNA melting, then it is no longer so clear to which of these lengths L should be compared. 
In order to tackle this more complex case, Binder et al 24j derived a method which is based 



on the use of an irrelevant variable u and an expression of the form 

f si n 9 = N- d F(tN yt ,uN y ") (22) 

After several approximations and a little bit of algebra, these authors obtain cy in the form 
of Eq. (|T9|) with, however, 

2d 
2(3 + 1 - 

Finally, the lengths that diverge at DNA melting are £ and (y), the average separation 
between paired bases. One might wonder whether L should not be compared to (y) instead 
of £. In order to check this hypothesis, let us denote by A the characteristic exponent for 
(y) {(v) ^ A )- Remember that if the external field is proportional to y then (y) is the order 
parameter m, so that A is equal to /3, the critical exponent for m (m oc t 13 ). Let us next 
express the singular part of the free energy in the form 

fsing = N' d F {tN- x ' x ) . (24) 
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Differentiating Eq. (1241) twice with respect to t, one again obtains cy in the form of Eq. 
f|T9|) with, however 

2 , 
p = ---d 

} (25) 

Table I shows the values of p and a calculated from the characteristic exponents reported 
in 2jl and Eqs (1201) . (l2Til . (1231) and (1251) . as well as adjusted values. These latter ones were 
obtained by varying p and o by hand in order that the plots of Cy/N p as a function of 
tN a are superposed for an interval of values of N as large as possible. By setting t = in 
Eq. (fT9|) . one sees that the maximum of cy scales as N p . p was therefore adjusted in the 
neighbourhood of the slope of the plot of the maximum of cy as a function of N (bottom plot 
of Fig. 1). On the other hand, a was adjusted in the neighbourhood of \jv. Examination 
of Table I indicates that the values of p and a obtained from Eqs. ( 120]) and ( |25i) compare 
well with the adjusted ones, while this is certainly not the case for the values obtained from 
Eqs. ( ]2Tj) and (123|) . Figs. 5 and 6 further show plots of Cy/N p as a function of tN a for, 
respectively, the DPB and JB models, and values of p and o obtained from Eq. (1211) (top 
plots) and adjusted ones (bottom plots). It is seen in the top plots that the curves with 
different values of N are far from being superposed for the values of p and a obtained from 
Eq. ( 12T1) . and the situation is still worse with Eq. ( 1231) . In contrast, the various curves 
are fairly well superposed for the adjusted values of p and a (see bottom plots of Figs 5 
and 6), as well as those obtained with Eqs. (120]) and fl25|) . An exception occurs for the 
curves corresponding to the largest values of N in the JB model (bottom plot of Fig. 6). 
Remember that the corresponding points also depart from the power law in the bottom 
plot of Fig. 1. The reason for this is that the TI method fails to give correct values of 
thermodynamical observables too close to the phase transition discontinuity. Examination 
of the top plots of Fig. 3 shows that sequences of length N = 10000 are still rather far 
from the thermodynamic limit for the DPB model, so that one needs not to worry about 
the effect of the discontinuity on TI calculations. In contrast, the bottom plot of Fig. 3 
and the two plots of Fig. 4 indicate that sequences of length N = 10000 have reached the 
thermodynamic limit for the JB model, so that the perturbative effect of the discontinuity 
becomes noticeable in TI calculations. 
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The fact that Eq. (I20I) leads to a correct superposition of the curves for different values 
of the sequence length N is the proof that the basic hypothesis of finite-size scaling theory 
is satisfied. Since, however, Eq. (1251) also leads to a correct superposition of the curves, it 
is, as expected, no longer clear to which diverging length (£ or (y)) L should be compared, 
both possibilities leading to a reasonable result. On the other hand, the fact that curves 
with different N are no longer superposed when Eq. (I2T]) is used to calculate p and a simply 
reflects the fact that Josephson's identity, 2 — a = ud, is not valid for these two models of 
DNA melting, a conclusion which was already arrived at in our preceding work 21] . Finally, 
the fact that curves also do not superpose when Eq. ( 1231) is used indicates that one of the 
several hypotheses made by the authors of Ref. 24j to arrive to these expressions is not 
satisfied for the DNA models, although it is not an easy task to tell which one(s) is(are) 
invalidated. Alternatively, Eq. (1231) can be straightforwardly derived from Eq. fl20l) by using 
Rushbrooke identity (a + 2(3 + 7 = 2) as well as Josephson's one. This latter identity being 
not valid, it comes as no surprise that Eq. ( |23l) leads to as poor a result as Eq. ( |2TI) . 



V. CONCLUSION 



To summarize, we modified the Transfer Integral method to adapt it to the calculation of 
thermodynamic quantities of finite sequences with open boundary conditions. Non-extensive 
quantities, like the average separation of paired bases (y) and the correlation length £, turned 
out to be the most tricky ones to evaluate. We then applied this modified procedure to the 
DPB and JB dynamical models, in order to clarify how the finiteness of the length of 
the sequence affects the phase transition that takes place at DNA melting temperature. 
We showed that the rounding of the transition that occurs when the size of the sequence 
decreases is clearly reflected in the temperature evolution of most quantities, including the 
specific heat cy, the correlation length £, and the average separation of paired bases (y). 
We next performed a finite-size scaling analysis of the two systems and showed that the 
singular part of the free energy can indeed be expressed in terms of an homogeneous function. 
However, since both £ and (y) diverge at the melting transition, it is no longer clear whether 
the argument of the homogeneous function should be (a power of) L/£ or Lj (y). Moreover, 
Josephson's identity is satisfied for none of the investigated systems, so that the derivation 
of the characteristic exponents p and a, which appear in the asymptotic expression of the 
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specific heat cy, requires some care. 

The Transfer Integral (TI) method appears as the only efficient numerical tool to study 
the thermodynamics of DNA melting in detail. In the formulation used here, it however 
applies only to homogeneous chains, while it is well established that the heterogeneity of 
real DNA molecules may smear out the discontinuity of the melting transition, just like 
the finiteness of the sequence does. Our next goal is therefore to overcome the technical 
difficulty associated with the application of the TI method to inhomogeneous chains and 
investigate the effect of heterogeneities on the phase transition at DNA melting. 

APPENDIX A: COMPUTATION OF THE CORRELATION LENGTH 

The static form factor is defined as 



S(q) 



and the correlation length is given by 



N 



^2(Vn ~ (Vn))e 



iqan 



n=l 



(Al) 



1 d 2 S{q) 



(A2) 



9=0 



2S(q) dq 2 

We stress that the statistical weight of the bases at the extremities n = 1, N is different from 
that of the other ones, so that they must be treated separately. We first write Eq. (lAlj) in 
the explicit form 



N N 

S ^ = E E <*jwW eiqa(n ~ m) ( A3 ) 

n=l m=l 

where 5y n = y n — (y n ). By isolating averages concerning extremity values, we get 



S(q) = (6y 2 ) + (6y 2 N ) + 2 {8y x 8y N ) cos[ga(iV - 1)] 
+ S ie iqa + S{e~ iqa + S N e lNqa + S* N e- lNqa + S mid 



(A4) 



where 
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N-l 

-iqam 



s 1 = Y<^ 5 y™) e ~ iq 

m=2 
N-l 

s N = Y,^y^) e ~ iqam ( A5 ) 



m=2 

N-l N-l 
m=2 n=2 



Defining 



dye -PVM(y)/2 y 2^ {y) 



we obtain the relations : 



(i)v(D J Ar^f^Af- if m > n 



(A6) 



(A7) 



Z t* [ A^ 1 A J n - m Af- n if m < n . 

According to the relations in Eqs. ( Till and (1A7I) . the summations in Eq. (1A3I) are just 
geometric series, which we evaluated formally in order to increase the speed of numerical 
calculations by a factor A^ 2 . After some tedious algebra, one obtains : 

S(q) = (Syl) + (5y%)+2 (5 yi 5y N ) cos[ga(iV-l)]+^ H i3 f f \ N) 

ij 1 Ji ^ l > 

+ S cosh^)- cos(ga) {Pv/wU) + C tJ f tJ (N)] cos[(iV - 2)qa] - P,jj)j{ \ - 1) - C<,-/«(2) 
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ijk 



-[DM2) + CMN - 1)] cos[(iV - l)qa] + [D^f^N) + C^-(l)] cos[qa}} 

2Mijk 

(1 - g(l, 1))(1 + ^(2) - 2^(1) cos(ga))(l + f jk (2) - 2f jk (l) cos(ga)) 



(A8) 



x {-0(3, 3) - 0(4, 4) - 0(2, 4) + g(N, N) + g(N, N + 2) + g(N + 1, N + 1) 

+ [0(2,3)-0(AT,Ar+l) + 20(3,4)-0(^+l,iV+2)+0(4,3)-0(iV+l,iV)-0(iV,iV+l)]cos(ga) 

+ [0(iV+l, iV+l)-(/(3, 3)] cos(2ga) + [0(2, N+l)-g(3, N+2)} cos[(A-3)ga] 

+ [0(4, N + 2)- 0(2, iV)] cos[(iV - 2)qa] + [g(3, N) - g(4, N + 1)] cos[(iV - l)qa}} 

2 



e 2(iqa-aij) _ e N(iqa-a i:j ) 



where 



\ giqa—aij 



Oii 



M,. 



Hi 



T- 



fij( n ) 
g(n,m) 



- In(r^) 



Cij = Tij - (01) Mij 



-notij —motjk 



(A9) 



Finally, S(q) must be derivated twice with respect to the wave vector q in order to get 
the correlation length from Eq. (IA2I) . 
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TABLE CAPTION 



Table I : Values of the coefficients p and a of Eq. ( fT9l) for the DPB and JB models, 
'irst four lines show the values calculated from the characteristic exponents reported in Ref. 



2l| and Eqs ( 1201) . ( I2T1) . ( |23i) and ( 1251) . The last line shows the values adjusted by hand in 



order that the plots of cy /N p as a function of tN a are superposed on an interval of values 
of N as large as possible (see bottom plots of Figs. 5 and 6). 
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FIGURE CAPTIONS 

Figure 1 : (color online) Log-log plots, as a function of the sequence length N, of 
the reduced critical temperature shift 1 — T c (N)/T c (oo) (top plot), the maximum of the 
correlation length £ (middle plot), and the maximum of cy (bottom plot), according to 
the DPB (squares) and JB (circles) models. £ is in units of the separation between two 
successive base pairs, and cy in units of k B . The solid and dash-dotted lines show the result 
of the adjustment of power laws against the calculated points. 

Figure 2 : (color online) Plot of the entropy per site s as a function of the rescaled 
temperature t(N) for an infinitely long chain (circles) and a sequence with N = 100 bp 
(squares), according to the DPB model (top plot) and the JB one (bottom plot), s is in 
units of k B . 

Figure 3 : (color online) Log-log plots of the specific heat per site cy as a function of 
the opposite —t(N) of the rescaled temperature for the DPB model (top plot) and the JB 
one (bottom plot) and seven values of the sequence length N ranging from 100 to oo. cy is 
in units of k B . Note that, at the thermodynamic limit of infinitely long chains, cy becomes 
infinite at the critical temperature but numerical limitations of the TI method prevent the 
observation of such divergence. 

Figure 4 : (color online) Log-log plots of the correlation length £ (bottom plot) and 
the order parameter (y) (top plot) as a function of the opposite —t(N) of the rescaled 
temperature for the JB model and seven values of the sequence length N ranging from 100 
to oo. £ is in units of the separation between two successive base pairs and (y) in units of 
the inverse 1/a of the Morse potential parameter. Although numerical limitations of the 
TI method prevent the observation of such divergences, £ becomes infinite at the critical 
temperature at the thermodynamic limit of infinitely long chains, while (y) becomes infinite 
at the critical temperature whatever the length N of the sequence. 

Figure 5 : (color online) Plots, for six values of the sequence length N ranging from 100 
to 10000, of Cy/N p as a function of tN a for the DPB model and values of p and o obtained 
from Eq. (12T1) (top plot) or adjusted by hand (bottom plot). 



19 



Figure 6 : (color online) Plots, for six values of the sequence length N ranging from 100 
to 10000, of Cy/N p as a function of tN a for the JB model and values of p and a obtained 
from Eq. (I2TI) (top plot) or adjusted by hand (bottom plot). 
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